|astro-ph/9806185 
September 10, 1999 
to appear in ApJ, 531, #2 (March 10, 2000) 



Supernova pencil beam survey 

Yun WangQ 

Princeton University Observatory 
Peyton Hall, Princeton, NJ 08544 
email: ywang@astro.princeton. edu 

ABSTRACT 

Type la supernovae (SNe la) can be calibrated to be good standard candles 
at cosmological distances. We propose a supernova pencil beam survey that 
could yield between dozens to hundreds of SNe la in redshift bins of 0.1 up 
to z = 1.5, which would compliment space based SN searches, and enable the 
proper consideration of the systematic uncertainties of SNe la as standard 
candles, in particular, luminosity evolution and gravitational lensing. We 
simulate SNe la luminosities by adding weak lensing noise (using empirical 
fitting formulae) and scatter in SN la absolute magnitudes to standard candles 
placed at random redshifts. We show that flux-averaging is powerful in reducing 
the combined noise due to gravitational lensing and scatter in SN la absolute 
magnitudes. The SN number count is not sensitive to matter distribution in the 
universe; it can be used to test models of cosmology or to measure the SN rate. 
The SN pencil beam survey can yield a wealth of data which should enable 
accurate determination of the cosmological parameters and the SN rate, and 
provide valuable information on the formation and evolution of galaxies. 

The SN pencil beam survey can be accomplished on a dedicated 4 meter 
telescope with a square degree field of view. This telescope can be used to 
conduct other important observational projects compatible with the SN pencil 
beam survey, such as QSOs, Kuiper belt objects, and in particular, weak lensing 
measurements of field galaxies, and the search for gamma-ray burst afterglows. 



1. Introduction 



Toward the end of the millennium, cosmology has matured into a phenomenological 
science. Observational data now dominates aesthetics in the evaluation of cosmological 
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models. Of fundamental importance is the determination of cosmological parameters, in 
particular, the Hubble constant H , the matter density faction Q m , and the density fraction 
contributed by the cosmological constant Observation of distant Type la supernovae 
(SNe la) has become an increasingly powerful means of measuring cosmological parameters 
( [Perlmutter et al. 1997| , 1998; [Riess et al. "199&1 ; |Schmidt et al. 1998Q , because SNe la 



can be calibrated to be good standard candles at cosmological distances. ([Riess, Press, & 
Kirshner 1995|) 



A type la SN is the thermonuclear explosion of a carbon-oxygen white dwarf in a binary 
when the rate of the mass transfer from the companion star is high. The SN explosion 
blows the white dwarf completely apart. The radioactive decay of the isotopes 56 Ni and 
56 Co is responsible for much of the light emitted. The SN lightcurve reaches a maximum 
about 15 days after the explosion and then declines slowly over years. A SN can outshine 
the galaxy in which it lies. 



Two independent groups ( [Riess et al. 1998| Perlmutter et al. 1999|) have made 



systematic searches for SNe la for the purpose of measuring cosmological parameters. 
Their preliminary results seem to indicate a low matter density universe, possibly with 
a sizable cosmological constant. Even though the uncertainty in these results is large, 
they clearly demonstrate that the observation of SNe la can potentially become a reliable 
probe of cosmology. However, there are important systematic uncertainties of SNe la as 
standard candles, in particular, luminosity evolution and gravitational lensing. To constrain 
the evolution of SN la peak absolute luminosities, we need a large number of SNe la at 
significantly different redshifts (low z and z > 1), which is not available at present. Both 
groups have assumed a smooth universe in their data analysis, although they include lensing 
in their error budgets. Since we live in a clumpy universe, the effect of gravitational lensing 
must be taken into account adequately for the proper interpretation of SN data. At present, 
the small number of observed high z SNe la prevents adequate modeling of gravitational 
lensing effects. 

In this paper, we propose a pencil beam survey of SNe la that could yield between 
dozens to hundreds of SNe la in redshift bins of 0.1 up to z = 1.5 (see §3) which would 
allow the proper modeling of gravitational lensing, as well as a quantitative understanding 
of luminosity evolution. Such a survey would yield a wealth of data which can be used 
to make accurate measurement of cosmological parameters and the SN rate, and provide 
powerful constraints on various aspects of the cosmological model. In §2, we consider the 
weak lensing of SNe la. We simulate SN la luminosities by adding weak lensing noise and 
scatter in SN la absolute magnitudes to standard candles placed at random redshifts. We 
show how flux- averaging reduces the combined noise of gravitational lensing and scatter in 
SN la absolute magnitudes. In §3, we show that the number count of SNe from a pencil 
beam survey is not sensitive to matter distribution in the universe; it can be used as a test 
of models of cosmology and SN progenitors, or to measure the SN rate accurately. In §4, we 
discuss the observational feasibility of the SN pencil beam survey. §5 contains conclusions. 
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2. Weak lensing of supernovae 

In a SN la Hubble diagram, one must use distance-redshift relations to make theoretical 
predictions. Unlike angular separations and flux densities, distances are not directly 
measurable, but they are indispensable theoretical intermediaries. The distance-redshift 
relations depend on the distribution of matter in the universe. 

In a smooth Friedmann-Robertson- Walker (FRW) universe, the metric is given by 
ds 2 = dt 2 — a 2 (t)[dr 2 /{l — kr 2 ) + r 2 (d6 2 + sin 2 6* <i0 2 )], where a(t) is the cosmic scale factor, 
and k is the global curvature parameter (Q^ = 1 — Q rn — = — k/H^). The comoving 
distance r is given by 



where "sinn" is defined as sinh if > 0, and sin if Qj, < 0. If — 0, the sinn and f^'s 
disappear from Eq. ([]]), leaving only the integral. The angular diameter distance is given by 
(1a{z) — r (^)/(l + z )i an d the luminosity distance is given by di(z) = (1 + z) 2 dA{z). 

However, our universe is clumpy rather than smooth. According to the focusing 
theorem in gravitational lens theory, if there is any shear or matter along a beam connecting 
a source to an observer, the angular diameter distance of the source from the observer is 
smaller than that which would occur if the source were seen through an empty, shear-free 
cone, provided the affine parameter distance (defined such that its element equals the 
proper distance element at the observer) is the same and the beam has not gone through 
a caustic. An increase of shear or matter density along the beam decreases the angular 
diameter distance and consequently increases the observable flux for given z. (Schneider, 
Ehlers, & Falco 1992) 

The observation of SNe la at z > 1 is important for the determination of cosmological 
parameters, but the dispersion in SN la luminosities due to gravitational lensing can become 
comparable to the intrinsic dispersion of SNe la absolute magnitudes because the optical 
depth for gravitational lensing increases with redshift. 



2.1. Direction dependent smoothness parameter 

If only a fraction a (known as the smoothness parameter) of the matter density is 
smoothly distributed, the largest possible distance (for given redshift) for light bundles 
which have not passed through a caustic is given by the appropriate solution to the following 
equation: 

+ ~an m {l + z) 5 D A = 0, (2) 



dD A 

dz 
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where g{z) = (1 + z)Vl + fi m ^ + ft A [(l + ^)~ 2 - 1]- ( |Kantowski T998] ) The Q A = form 
of Eq. (§) has been known as the Dyer-Roeder equation, ( pyer fc Roeder 1973|, |Schneider et 
al. 1992| ) 



Fig. 1(a) shows magnitude versus redshift for the three cosmo logical models considered 
by Riess et al. (1998), SCDM (Q m = 1, Q A = 0), OCDM (Q m = 0.2, fi A = 0), and 
ACDM (fi m = 0.2, £7 a = 0.8). For each cosmological model, the upper curve represents 
the completely clumpy universe (empty beam, 5 = 0), while the lower curve represents the 
completely smooth universe (filled beam, a = 1). Fig. 1(b) shows the same models relative 
to smooth OCDM ( filled beam, 5 = 1), the middle curve for each model now represents a 
universe with half of the matter smoothly distributed (half-filled beam, a = 0.5). Clearly, 
at z > 1, there is degeneracy of distances in a flat clumpy universe and an open smooth 
universe, and also in an open clumpy universe and a flat smooth universe with a sizable 
cosmological constant, as has been noted by a number of previous authors. ( |Kantowskl 



1998 




Linder 1998, 


Holz 1998) 



We can generalize the angular diameter distance D^(z) by allowing the smoothness 
parameter a to be direction dependent, i.e., a property of the beam connecting the observer 
and the standard candle. The smoothness parameter a essentially represents the amount 
of matter that causes weak lensing of a given source. Since matter distribution in our 
universe is inhomogeneous, we can think of our universe as a mosaic of cones centered on 
the observer, each with a different value of a. This reinterpretation of a implies that we 
have a > 1 in regions of the universe in which there are above average amounts of matter 
which can cause magnification of a source. ( |Wang 1999a| ) 



In order to derive a unique mapping between the distribution in distances and the 
distribution in the direction dependent smoothness parameter for given redshift z, we define 
the direction dependent smoothness parameter a to be the solution of Eq.(Q) for given 
distance D^(z). 

At given redshift z, the magnification of a source can be expressed in terms of the 
apparent brightness of the source f(a\z), or in terms of the angular diameter distance to 
the source DA{oi\z): 

Da{ch = l\z) 



f(a\z) 



D A (a\z) 



(3) 



f(a = l\z) 

where f(a = l\z) and D^a = l\z) are the flux of the source and angular diameter distance 
to the source in a completely smooth universe (filled beam), and a is the direction dependent 
smoothness parameter. Since distances are not directly measurable, we should interpret 
Eq.(|3[) as defining a unique mapping between the magnification of a standard candle at 
redshift z and the direction dependent smoothness parameter a at z; a parametrizes the 
direction dependent matter distribution in a well-defined manner. 



From the magnification distributions of standard candles at various redshifts, p(/j,\z), 
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with z =0.5, 1, 1.5, 2, 2.5, 3, 5, found numerically by Wambsganss et al. for Vt m = 0.4, 

= 0.6 ([Wambsganss et al. 1997| , [Wambsganss 1999|) , Wang (1999a) has obtained simple 
empirical fitting formulae for the distribution of a: 

21 



p(a\z) = C n 



exp 



w a q 



(4) 



where C norm , a pea k, w, and q depend on z and are independent of a. They are given by 
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(5) 



C n0 rm{ z ) is the normalization constant for given z. The parameter a pea k{z) indicates 
the average smoothness of the universe at redshift z, it increases with z and approaches 
& P eak(z) = 1 (filled beam) at z = 5; the parameter w(z) indicates the width of the 
distribution in the direction dependent smoothness parameter a, it decreases with z. The z 
dependences of a pea k{. z ) and w(z) are as expected because as we look back to earlier times, 
lines of sight become more filled in with matter, and the universe becomes smoother on 
the average. The parameter q(z) indicates the deviation of p(d\z) from Gaussianity (which 
corresponds to q — 0). 

Models with different cosmological parameters should lead to somewhat different 
matter distributions p(a\z). In the context of weak lensing of standard candles, we expect 
the cosmological parameter dependence to enter primarily through the magnification /i to 
direction dependent smoothness parameter a mapping at given z (the same a corresponds 
to very different fi in different cosmologies). 



2.2. Flux averaging of SN luminosities 

Gravitational lensing noise in the Hubble diagram can be reduced by appropriate flux 
averaging of SNe la in each redshift bin. Because of flux conservation, the average flux of 
a sufficient number of SNe la at the same z from the same field should be the same as the 
true flux of the SNe la without gravitational lensing if the sample is complete. 

It is convenient to compare the distance modulus of SNe la, /io, with the theoretical 
prediction 
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where di{z) is the luminosity distance. Before flux-averaging, we convert the distance 
modulus fio(zi) of SNe la into "fluxes", f(zi) = lO - ^ 2 ^/ 2 - 5 . We then obtain "absolute 
luminosities", {C(zi}}, by removing the redshift dependence of the "fluxes", i.e., 



£(zj) = 4:7id 2 L (zi\H ,n m ,n A ) f(zi 



(7) 



where (H , Q m , are the best-fit cosmological parameters derived from the unbinned data 
set {f(zi)}. We then flux-average over the "absolute luminosities" in each redshift bin. 
The set of best-fit cosmological parameters derived from the binned data is applied to the 
unbinned data {f(zi}} to obtain a new set of "absolute luminosities" {£«}, which is then 
flux-averaged in each redshift bin, and the new binned data is used to derive a new set of 
best-fit cosmological parameters. This procedure is repeated until convergence is achieved. 
This iteration should lead to the optimal removal of gravitational lensing noise and the 
accurate determinations of the cosmological parameters. Wang (1999b) has applied this 
method to analyze the combined data from the two groups ( [Riess et al. 1998] , |Perlmutter 
et al. 1999| ). 



To illustrate how flux- averaging can reduce the dispersion in SN la luminosities caused 
by weak lensing, let us simulate the data by drawing Nsn random points from the redshift 
interval [^i,^]) each point represents a standard candle. We add weak lensing noise by 
giving each standard candle a direction dependent smoothness parameter a (corresponding 



to fi = [Da((x = 1)/Da{c()] 2 ) drawn at random from the distribution p(a\z) given in 
The scatter in SN la absolute magnitudes, Am abm) can be written as 



52.1. 



Am, 



abm 



Am int + Am 



obs ) 



where Amj„ ( is the intrinsic scatter and Am^ is the observational noise. We assume that 
both intrinsic scatter and observational noise are Gaussian distributed in magnitude, with 
dispersions ai nt and a b s respectively. Then the total scatter in SN la absolute magnitudes 



is also Gaussian distributed, i.e., 



p(Am abm ) 



1 



2t Cabm 



exp 



(Am abm f 

^ a abm 



(9) 



with a abm = yafnt + cr^ bs . We take 
extracted from the data is 



0.20. The absolute luminosity of each SN la 



C(z) = 10- Am °»^ 2 - 5 Li(a\z)C int (z) 7 

= iO" Am ^/ 2 - 5 fj,(a\z) {lO~ Am ™ t/2 - 5 £(a = l\z)} 
= C(a = l\z) n(a\z) iO" Am ^/ 2 - 5 . 



(10) 



We have used C int (z) = i(r A "W 2 - 5 c{a = \\z). 
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For each SN la, the total noise is 



Let us average the fluxes of all SNe la in the redshift bin [z®, z®]' 



The flux averaged noise is 



1 Nsn 

' £ £(*)• (12) 



(AmU = "2-5 log (zja^) • («) 

where 

\^N SN _ 

-=%^. (14) 

™ SN 

Table 1 lists the means and dispersions (in the form of meanidispersion) of Am 
(which are (Am) and a = y ([Am — (Am)] )), and {Am) avg (which are ((Am) avg ) and 

a avg = \J ([(Am) avg — ((Am) aV g)] 2 ) ) for various redshift bins with N S n=2, 4, and 9 SNe la 
in each bin, for 10 4 random samples. We have taken Vl m = 0.4, fi A = 0.6. 



Table 1: Means and dispersions of Am and {Am) avg for a ACDM model. 



(Am) ± a 



N SN =2 



N SN =4 



N SN =9 



[0.5,0.6] 
[1,1.1] 

[1.5,1.6] 
[1.5,2] 
[2,2.5] 
[2.5,3] 
[3,3.5] 
[3-5,4] 
[4,4.5] 
[4.5,5] 



0.001 ±0.200 
0.002 ±0.204 
0.003±0.210 
0.003±0.213 
0.005 ±0.218 
0.006 ±0.223 
0.007 ±0.228 
0.007 ±0.233 
0.008 ±0.238 
0.008±0.243 



-0.007±0.141 
-0.007 ± 0.145 
-0.006± 0.149 
-0.006± 0.151 
-0.005 ±0.155 
-0.005±0.159 
-0.004± 0.162 
-0.004±0.166 
-0.004±0.170 
-0.004± 0.174 



-0.012 ±0.101 
-0.013 ± 0.103 
-0.012 ± 0.107 
-0.012 ± 0.108 
-0.012 ± 0.111 
-0.012±0.114 
-0.012±0.117 
-0.012±0.119 
-0.012±0.122 
-0.012± 0.125 



-0.015±0.067 
-0.015 ± 0.069 
-0.015 ± 0.071 
-0.015 ± 0.072 
-0.015 ± 0.074 
-0.015±0.076 
-0.015± 0.078 
-0.015 ± 0.080 
-0.015± 0.082 
-0.016± 0.084 



The dispersion decreases roughly as \/^/N S n- The dispersion would reduce by 30% 
if the sample contains 2 SNe la; 50% if the sample contains 4 SNe la. Even though 
gravitational lensing noise increases with redshift, the combined gravitational lensing and 
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SN la absolute magnitude scatter noise in the redshift interval z = [4.5, 5] can be reduced 
to the same level as in the absence of lensing by flux averaging over two SNe la. 

Note that the flux averaged luminosities are biased towards slightly higher luminosities. 
This is as expected, because we have assumed that the intrinsic scatter and observational 
noise in the SN la absolute magnitude are Gaussian in magnitude. It is straightforward 
to show that the mean of lcr Amabm / 2 - 5 is exp [(In 10/2. 5) 2 c 2 bm /2], which corresponds to 
a bias of — o\ hm In 10/5 ~ —0.018 for a a bm = 0.2. If we assume that the intrinsic scatter 
and observational noise are Gaussian in luminosity, the flux averaged luminosities become 
unbiased. 



3. Supernova number count 



The SN number count is most sensitive to the SN rate as a function of z, which depends 
on the specific rate of SNe, as well as the number density of galaxies and their luminosity 
distribution. The frequency of SNe is a key parameter for describing the formation and the 
evolution of galaxies; the winds driven by SNe tune the energetics, and their production 
of metals determines the chemical evolution of galaxies and of clusters of galaxies (|Ferrini| 
fe Poggianti 1993| , |Kenzini et al. 1993| ). The SN II rate is related (for a given initial 
mass function) to the instantaneous stellar birthrate of massive stars because SNe II have 
short-lived progenitors; the SN la rate follow a slower evolutionary track, and can be used 
to probe the past history of star formation in galaxies. Accurate measurements of the SN 
rates at intermediate redshifts are important for understanding galaxy evolution, cosmic 
star formation rate and the nature of SN la progenitors. ( [Madau 1998| , |Ruiz-Lapuente & 
Canal 1998] , |Sadat et al. 1998| , [Yungelson fc Livio 1998] [Madau et al. 1998|) 



The SN rates are very uncertain at present due to the small number of SNe discovered 
in systematic searches. (|van den Bergh fc McClure 1994] , Pain et al. 1996|) . Kolatt 
& Bartelmann (1997) have estimated the SN la average rate per proper time unit per 
comoving volume to be 



n SN (z) = [A + Bz] (lOOyrT^/i^Mpc)- 



A = 0.0136,5 = 0.067, 



(15) 



for q = 0.5. Changing q to lower values leads to lower comoving densities but higher 
luminosities of galaxies at high z ([Lilly et al. 1995|) . Eq.(|15D has been derived assuming 
that all SNe reside in galaxies, neglecting the redshift dependence of the specific SN rate 
(number per unit luminosity per unit time), and using the number density of galaxies 
(as function of redshift) and the Schechter function parameters derived from the Canada 
France redshift survey ( |Lilly et al. 1995] ), the APM survey ( [Loveday et al. 1992| ), and the 
AUTOFIB survey ( [Ellis et al. 1996|) . We use Eq.(^) for all cosmological models considered 
in this paper, because it is a very crude and conservative estimate, and we use it for the 
purpose of illustration only. 
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The expected number of SNe la in a field of angular area 8 2 for an effective observation 
duration of At up to redshift z is 



dr p (z') 2 n SN (z') 
r \z ) 




(16) 



0.0136 + 0.067^' 



N(z) 



(i + z yn m {i + z'f + Ov + + z'f 

where dr p = —cdt is the proper distance interval, r is comoving distance, and the factor 
(1 + z)^ 1 accounts for the cosmological time dilation. Note that the number counts 
fundamentally probe a different aspect of the global geometry of the universe than do the 
distance measures flCarroll, Press, fc Turner 1992| ). 

Fig.2 shows the number of SNe la expected per 0.1 redshift interval as function of 
redshift for the same three cosmological models as in Fig.l, SCDM (Q m = 1, Q\ = 0), 
OCDM (Q m = 0.2, Q A = 0), and ACDM (£l m = 0.2, fi A = 0.8). For a one square degree 
field, and an effective observation duration of one year, the total numbers of expected SNe 
la are 464, 899, and 1705 for SCDM, OCDM, and ACDM respectively for z up to 1.5. 

The number of SNe which are strongly lensed by galaxies is given by 

Nlensed{z) = [ dz'r(z')^-, 

Jo dz' 



3.81 x 10 




0.0136 + 0.067,2' 



1 + z 



n k (i + z 



'\2 



(17) 



We have used the optical depth for gravitational lensing by galaxies ( [lurner 1990| , [Kukugita 
fe Turner 1991 



t(z) 



F 
30 



[l + z)d A (z) 
cH, 1 



'11 



where F parametrizes the gravitational lensing effectiveness of galaxies (as singular 
isothermal spheres). Fig. 3 shows the number of strongly lensed SNe as function of survey 
depth z, for the same cosmological models as in Fig.l and Fig.2. For a one square degree 
field, and an effective observation duration of one year, the total numbers of strongly lensed 
SNe la are 0.2, 0.6, and 1.8 for SCDM, OCDM, and ACDM respectively for z up to 1.5. 

Fig.2 and Fig. 3 show that the SN number counts from a pencil beam survey can be 
used to measure the SN rate at high redshifts and perhaps to probe cosmology. Because 
of the large number of SNe in each redshift bin and the smallness of gravitational lensing 
optical depth, the SN number count should be insensitive to matter distribution in the 
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universe, and should therefore provide a robust probe of the cosmological model. Note 
that the strongly lensed SNe can be easily removed from the survey sample, because they 
appear as unusually bright SNe. The number of strongly lensed SNe is very sensitive to the 
cosmological model and might be used to further constrain the cosmological model. But 
given the small number of strongly lensed SNe la expected in any realistic observational 
program, their usefulness may be limited. 

The SN number counts provide a combined measure of the cosmological parameters 
and the SN rate. Fig. 4 shows the parameter dependence of the SN la number count per 
0.1 redshift interval (z — 0.05, z + 0.05) as function of redshift z. Note that the dependences 
of the SN number count on the SN rate (parametrized by A and B) and the cosmological 
parameters are not degenerate and can be differentiated in principle. In practice, we are 
ignorant of the functional form of the SN rate as a function of z. Hence, we may apply 
the measurements of the cosmological parameters as priors to the SN number counts to 
obtain a direct measure of the SN rate in the universe for z up to 1.5, which can be used 
as a powerful constraint on the cosmological model. For a one square degree field, and 
an effective observation duration of one year, the SN rate per 0.1 redshift interval can 
be determined to 14-16%, 9-11%, 7-8% at 1 < z < 1.5 for SCDM, OCDM, and ACDM 
respectively. 



4. Observational feasibility 

A pencil beam survey would be efficient in the discovery of SNe at z > 1 through the 
combination of data from successive nights and the comparison of the latest frame of images 
with all previous frames. 

SNe la are quite faint at z ~ 1.5, with AB magnitude in the I band of Iab ~ 26. 
Because of the UV suppression due to line blanketing and the apparent IR suppression 
in the rest frame SN la spectrum, one should use a passband that corresponds to the 
wavelength range of 3000A — lOOOOA in the SN rest frame; this means using the I, J, H, 
or K passbands to observe the z > 1 SNe. SN searches from space (HST and NGST) are 
limited by small fields of view, a large scale SN search is only possible from the ground, 
where one is limited to the I band by the atmosphere. Here we discuss the observational 
feasibility of a pencil beam survey of SNe la up to z ~ 1.5 in terms of exposure times in the 
I band; although multiple band photometry should be obtained to constrain extinction and 
evolution of the SNe. 

Using values for the photometric parameters from the Sloan Digital Sky Survey (1 
arcsec seeing, effective sky brightness of 20.3 mag/arcsec 2 in the I band, seeing-dominated 
PSF, etc), we find that the exposure time for a point source with AB magnitude of Iab can 
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be written as 

t = 13.94hours (^f lO - 8 ^*' 26 ) (19) 

where S/N is the signal-to-noise ratio, D is the aperture of the telescope. The above 
equation shows that the supernova pencil beam survey can be accomplished from a modest 
dedicated 4 meter telescope, which can image the same fields (two pencil beams should be 
observed to keep the fields close to zenith) every night, which can lead to the discovery of 
SNe la up to z — 1.5 via appropriate combination of data from successive nights, and the 
light curves of SNe la at z < 1.5. It should be feasible to monitor the faintest SNe la from 
the Keck 10 meter telescope or the HST. 

Spectra of SNe are required to determine whether they are type la. SNe la have a 
Si II absorption line at ~ 4000 A that may be used to identify the 1 < z < 1.5 SNe la 
in ground-based observations. The follow-up spectroscopy can be attempted on the Keck 
Low- Resolution Imaging Spectrometer (LRIS). Assuming that we use the 300 grooves/mm 
grating which provide dispersions of 4.99 A per 48 /iin (2 pixels on the CCD), we find the 
exposure time for 0.5 arcsec seeing to be 

t = 28 hours ^Vf—^Vo a4[2( ^ fl - 26) - (/ -- 21)1 , 
\ 3 /VI arcsec / 

= 4.44 hours ( ( -^-] 10 MP(i^-25)-<i#-2i)] 
\ 3 /VI arcsec / 

= 12.33hours ( (— — ) 10 ^^)-(^-^, (20) 
V 5 J VI arcsec / 

where W is the width of the slit, and I^b is the sky brightness in the I band in units of 
mag/arcsec 2 . Note that the Si II absorption feature at ~ 4000 A is not deep enough to be 
useful in a very noisy (S/N = 3) spectrum. Clearly, spectroscopic follow-up of the z ~ 1.5 
(Iab ~ 26) SNe la will require substantial observational resources; NICMOS or NGST will 
be more suitable for the spectroscopy of the SNe la at the highest redshifts discovered 
from the ground. The Keck LRIS can be used to obtain the spectra of the SNe la at more 
modest redshifts (z ~ 1.3, Iab ~ 25). 

A dedicated 4 meter telescope with a square degree field of view can be used to conduct 
other important scientific projects compatible with the SN pencil beam survey, such as 
QSOs, Kuiper belt objects, and in particular, weak lensing and the search of gamma-ray 
burst (GRB) afterglows. 

Weak lensing is a powerful tool in mapping the mass distribution in the universe. The 
large field of view and the depth of a pencil beam survey would be ideal for weak lensing 
measurements of field galaxies, which can be used to constrain the large scale structure in 
the universe. 
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GRBs are perhaps the most energetic astrophysical events in the universe. Currently, 
there are many competing theories to explain GRBs. GRB afterglows contain valuable 
information, a statistically significant sample of GRB afterglows can provide strong 
constraints on the GRB theory. If beaming is involved in GRBs, we expect most of the GRB 
afterglows not to be associated with observable bursts. Schmidt et al. (1998) have observed 
a few optical transients which were too short in duration to be SNe; these could be GRB 
afterglows. Since present observational data seem to indicate that GRBs have associated 
host galaxies, the detection of host galaxies associated with short optical transients would 
support the interpretation of the latter as GRB afterglows. Since the GRB host galaxies 
are typically fainter than R = 25, a pencil beam survey would be ideal in detecting these 
host galaxies of candidate GRB afterglows. 



5. Summary and discussion 

We have proposed a pencil beam survey of SNe la which can yield from tens to 
hundreds of SNe la per 0.1 redshift interval for z up to 1.5, which would enable the 
quantitative consideration of the systematic uncertainties of SNe la as standard candles, 
in particular, luminosity evolution and gravitational lensing. Using the Perlmutter et al. 
"batch" search technique repetitively over the same field in the sky, the pencil beam survey 
would be efficient in the discovery of SNe at z > 1 by allowing the comparison of the latest 
frame of images (which may consist of combined data from successive nights) with all 
previous frames. 

The direct product of such a survey is the SN number count as a function of redshift 
(see §3), which is a combined measure of the cosmological parameters and the SN rate. 
When the measurements of the cosmological parameters are applied as priors to the number 
count, we obtain a direct measure of the SN rate, which is a key parameter in the formation 
and evolution of galaxies. The non-type-la SNe discovered by the pencil beam survey may 
be comparable to the type la SNe in number. 

The most important and straightforward application of the data from the SN pencil 
beam survey is to reduce the gravitational lensing noise in a SN la Hubble diagram via 
flux averaging. We have simulated SN la luminosities by adding weak lensing noise (using 
empirical fitting formulae given by Wang 1999a) and scatter in SN la absolute magnitudes to 
standard candles placed at random redshifts. We have shown that flux-averaging is powerful 
in reducing the combined noise of gravitational lensing and SN la absolute magnitude 
scatter (see §2). Because of the non-Gaussian nature of the luminosity distribution of SNe 
la at given z due to weak lensing, the large number of SNe la in a given redshift interval at 
high z is essential for the proper modeling and removal of the gravitational lensing effect. 
The SN la luminosity distribution in each redshift interval can be used to constrain the 
cosmological model (in particular, the fraction of matter in compact objects) by comparison 
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with predictions of numerical ray-shooting. ([Holz fc Wald 1998|) 

The completeness of the SN la sample determines the effectiveness of the removal of 
gravitational lensing noise from the SN la Hubble diagram and the amount of information 
contained in the SN la luminosity distribution in each redshift interval. Note that the 
magnitude limit of the survey can lead to observational bias against the most distant 
demagnified SNe, therefore, the SNe which are close to the magnitude limit of the survey 
should not be used to probe cosmology in the manner described in this paper. To ensure the 
maximum usefulness of the data, scrupulous attention will have to be paid to photometric 
calibration, uniform treatment of nearby and distant samples, and an effective way to deal 
with reddening. ( [Riess et al. 1998|) 



Note that although we can remove or reduce the effect of gravitational lensing on 
the SN la Hubble diagram, other systematics can affect the observed luminosity of SNe 
la. For example, grey dust, an evolution of the reddening law, or evolution in the peak 
absolute luminosity of SNe la. While it is possible to constrain grey dust and determine 
dust evolution through multi-band photometry at significantly different redshifts, the 
dimming with z in peak absolute luminosity of SNe la is degenerate with the effect of low 
matter density. QAguirre 1999 , [Riess et al. 1999 . Drell. Lorcdo, fc Wasscrman 1999 , Wang 



1999b|) Evolution will remain a caveat in the usage of SNe la as cosmological standard 



candles, unless one can somehow correct for the effect of evolution. It is critical to obtain 
up to hundreds of SNe la at z > 1 to constrain luminosity evolution, because we expect 
the luminosity evolution and low matter density to affect the distance modulus of SNe la 
through different functionals of z, which should become distinguishable at z > 1. 

We have proposed a survey cutoff of z = 1.5 mainly for two reasons. First, going to 
higher redshift makes obtaining the spectra of the SNe (which are needed to distinguish 
different types of SNe) practically impossible from the ground; even obtaining the spectra 
of z ~ 1.5 SNe la may prove impractical from the ground (as it would require several clear 
nights per spectrum on the Keck, see §4). Observers have already demonstrated that SNe la 
up to z = 1 can be found in the ongoing searches. flGoobar fc Perlmutter 1995| , |Garnavich 



et al. 199~S| , |Perlmutter et al. 1999| ) At z ~ 1.5, the SNe la should be ~ 2.5 magnitudes 



fainter than at z ~ 1, it should be possible to discover them on the ground through deep 
imaging, and the follow-up spectroscopy can be done on NICMOS or NGST (see §4). 
Second, predicted rest-frame SN la rate per comoving volume as function of redshift seems 
to peak at z ~ 1. ([Yungelson fc Livio 1998| , |Sadat et al. 19981) A pencil beam survey of 



SNe up to z — 1.5 will enable the accurate determination of the SN rate as function of 
redshift in the redshift region important for studying cosmic star formation rate and the 
SN la progenitor models. Although the Next Generation Space Telescope can detect SNe 
at as high redshifts as they exist QStockman et al. 1998|) , the estimated rate of detection is 



of order 20 SNe II per 4x4 arcmin 2 field per year in the interval 1 < z < 4 (Madau et al. 



1998|) , and the detection rate of SNe la is likely smaller. Thus a ground-based pencil beam 



survey of SNe is essential to complement the space-based SN searches. 
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The most challenging aspect of the SN pencil beam survey is obtaining spectra for the 
SNe la at redshifts close to 1.5. Instead of waiting for future space equipments, we may find 
innovative ways of obtaining spectra from the ground. The referee has pointed out that 
since we can find multiple SNe at the same time on the same one square degree field (the 
number of SNe depends on cosmology), it may be possible to get multiple spectra at once 
with fibers. 

The nominal numbers we have used for the proposed SN pencil beam survey, a one 
square degree field, a depth of z — 1.5, and an effective observation duration of one 
year (which is equivalent to several years of actual observation), are optimistic but not 
implausible. The large sky coverage and the long effective observation duration will 
probably require a large consortium of existing and new SN search teams through, for 
example, a dedicated 4 meter telescope which can be used at the same time for other 
important observational projects compatible with the SN pencil beam survey (see §4). The 
goal of going up to z = 1.5 in spectroscopy will require support from the Keck LRIS and 
NICMOS/NGST. We conclude by noting that a SN pencil beam survey can yield enormous 
scientific return. The observational efforts directed towards a SN pencil beam survey should 
be very rewarding. 
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Fig. 1. — (a) Magnitude versus redshift for three cosmological models, SCDM (Q m = 1, 
tt A = 0), OCDM (Q m = 0.2, tt A = 0), and ACDM (Q m = 0.2, fi A = 0.8). For each 
cosmological model, the upper curve represents the completely clumpy universe (empty 
beam, a — 0), while the lower curve represents the completely smooth universe (filled beam, 
a = 1). (b) The same models relative to smooth OCDM (filled beam, a — 1), the middle 
curve for each model represents a universe with half of the matter smoothly distributed 
(half-filled beam, a = 0.5). 

Fig. 2. — The number of SNe la expected per 0.1 redshift interval as function of redshift 
for three cosmological models, SCDM (fi m = 1, Q A = 0), OCDM (Q m = 0.2, fi A = 0), and 
ACDM (Q m = 0.2, f2 A = 0.8), for a one square degree field, and an effective observation 
duration of one year. 

Fig. 3. — The number of strongly lensed SNe as function of survey depth z, for the same 
cosmological models as in Fig.l and Fig. 2, for a one square degree field, and an effective 
observation duration of one year. 

Fig. 4. — The parameter dependence of the SN la number count per 0.1 redshift interval as 
function of redshift. Note that the dependence of the SN number count on the SN rate is 
parametrized by A and B. 
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Fig. 1. — (a) Magnitude versus redshift for three cosmological models, SCDM (Q m = 1, 
tt A = 0), OCDM (fi m = 0.2, tt A = 0), and ACDM (Q m = 0.2, fi A = 0.8). For each 
cosmological model, the upper curve represents the completely clumpy universe (empty 
beam, a — 0), while the lower curve represents the completely smooth universe (filled beam, 
a = 1). 
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Fig. 1. — (b) The same models relative to smooth OCDM (filled beam, a — 1), the middle 
curve for each model represents a universe with half of the matter smoothly distributed 
(half-filled beam, a = 0.5). 
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Fig. 2. — The number of SNe la expected per 0.1 redshift interval as function of redshift 
for three cosmological models, SCDM (fi m = 1, Q A = 0), OCDM (Q m = 0.2, Q A = 0), and 
ACDM (Q m = 0.2, f2 A = 0.8), for a one square degree field, and an effective observation 
duration of one year. 
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Fig. 3. — The number of strongly lensed SNe as function of survey depth z, for the same 
cosmological models as in Fig.l and Fig. 2, for a one square degree field, and an effective 
observation duration of one year. 
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Fig. 4. — The parameter dependence of the SN la number count per 0.1 redshift interval as 
function of redshift. Note that the dependence of the SN number count on the SN rate is 
parametrized by A and B. 



